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' We study the phase diagram of glassy systems in presence of an attractive coupling among real 

replicas. We find competition among a localized and a delocalized phase, that are separated by a 
coexistence line as in ordinary first order phase transitions. The coexistence line terminates in a 
critical point. We present numerical simulations for binary glasses in which show that this scenario 
is realized. 
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I. INTRODUCTION 



Supercooled liquids show a dramatic increase of the relaxation time as the temperature is lowered. The glassy 
transition is met at the temperature T g where flowing becomes not observable and the systems get out of equilibrium. 
One of the most suggestive ideas to rationalize this behavior can be found in the Gibbs-DiMarzio theory and its 
1 elaborations [Q. This relates the observed finite time singularity to a thermodynamic transition due to an entropy 
crisis, which would be observable only in the infinite time limit. Qualitatively the picture goes as follows. Still in 
the liquid phase, around a crossover temperature T c , the time scales for local motion of the molecules (vibrations) 
and large scale motion (diffusion) become widely separated, as if the system remained trapped for long times in local 
minima of the free energy. The entropy E associated to large scale motion, often called configurational entropy or 
complexity, can be related to the number of free-energy minima AT by the formula M = exp(iVE(T), (N is the number 
of particles). As a function of temperature, E is supposed to vanish at a finite temperature T s , which is the ideal point 
of thermodynamic glass transition. The relaxation time is then related to the escape rate from the free-energy minima 
and one can argue in favor of Vogel-Fulcher like relations for the relaxation time as a function of the temperature 
with a divergence at T s . As it was first realized by Kirkpatrick and Thirumalai Q , the Gibbs-DiMarzio scenario 
is exactly implemented in a large class of infinite range disordered models, with the difference that the times needed 
y—( • to escape from local equilibrium states (and the corresponding free-energy barriers) diverge when the volume of the 
system goes to infinity. The value of T c , which is to a large extent arbitrary in real systems, can be sharply defined 
G\ ■ in the mean-field limit. In fact this is the temperature where the mode coupling theory [Q, which is exact in these 
models, shows a divergent relaxation time. On the contrary in short range systems T c signals a change in behavior 
but we cannot assign to it any sharply defined value. 

It has been recently shown that if the models are generalized by introducing two coupled replicas of the same 
system one finds that T c corresponds to the edge of a metastability region j5j. In the same way the complexity is 
related to the difference of free-energy in the stable and in the metastable phase. In this paper we show how in the 
framework of coupled replicas the glass transition can be described as an ordinary phase transition. Enlarging the 
space of the parameters to include the coupling among replicas we find a first order transition line, terminating in a 
critical point. Although our analysis is based on mean-field theory, we will see that, as in ordinary first order phase 
• i-H , transition, the Maxwell construction will allow to extract the qualitative features of the phase diagram of real systems. 
A sketch of these results has appeared in ref. j6|. To submit to test our picture in realistic systems, we have simulated 
coupled replicas of binary mixtures with repulsive interaction |t],|^] . These are known to vitrify for some values of the 
parameters defining the model. The results of Monte Carlo simulations strongly support the theoretical picture. 

This paper is structured as follows. In section II we present some general considerations on the effect of coupling 
replicas and we predict the behavior of a glass in the presence of two coupled replicas of the same system. We have to 
distinguish two different case: the quenched and the annealed one, which have different properties. In section III we 
show that the previous arguments are indeed correct in a soluble model for the glassy transition, the p-spin spherical 
model. In section IV we present the results of the numerical simulations for the binary mixtures. Finally in the last 
section we present our conclusions. 
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II. COUPLING REPLICAS 



A. The quenched case 



In this section we describe the construction of an effective potential in systems where long range order is absent, 
but that can remain stuck for a long time in metastable states. Generically, one can expect that all the relevant 
metastable states at a given temperature are equivalent as far as their macroscopic characteristics are concerned. 
Consequently, one can not identify any intrinsic order parameter allowing to distinguish one state from the others. 
In this situation it is possible to use as order parameter a degree of similarity among different points in configuration 
space. The procedure is common in spin glass theory where the "overlaps" among different replicas appear as natural 
order parameters when one averages over the quenched disorder. Here we discuss an effective potential as a function 
of the overlap for generic systems, which may or may not, as structural glasses, contain quenched disorder. 

Let us describe the construction in the case of a systems composed by only one type of particles with coordinates 
Xi, for i = 1, JV; the generalization to many kind of particles is trivial. We consider two replicas of the same system, 
with coordinates x and y respectively, in an asymmetric relation. The replica y is a typical configuration distributed 
according to the Boltzmann-Gibbs law with the original Hamiltonian of the system H(y) at a temperature T', and 
does not feel any influence from the replica x. The replica x, instead, feels the influence of the replica y, and for fixed 
value of y, thermalizes at a temperature T with a Hamiltonian 

N 

H e (x\y)=H(x)-eJ2 w ( x i-yk) (1) 

z,fc=l 



The function w is different from zero only at short distance, an example is w(x) = 1 if \x\ < a and w(x) = 
if \x\ > a. An interesting behavior is found when the value of a is smaller that the typical interatomic distance 
(e.g. a = .3 atomic distances). The quantity q(x,y) = N^ 1 k=i n w ( x i ~ Vk) measures then the similarity of 
the two configurations, and would be close to one when the two replicas stay in similar configurations. Using the 
same terminology as in spin glasses J9| |ll[ q can be called the overlap of the two configurations. For positive e the 
x variables feel then a potential which pushes them near to the y variables. We can define a free-energy for the x 
variables in presence of the quenched configuration y as: 



F(T,e,y) = (Np) 



Mn / dxexp{-/3H(x) + fie ^ w ( x i 
y i,k=l 



Vk)} 



(2) 



This quantity should be self-averaging with respect to the distribution of the y and can therefore be computed as 

fdycM-P'H(y))F(T,e,y) 



F Q (T,T',e) = 



fdyeiq>(-{3>H(y)) 



(3) 



The temperature T' of the reference configuration y can be equal or different from that of the x configuration (T). 

The free-energy (^) is a well defined function that one can envisage to evaluate analytically or numerically. Its 
practical analytic evaluation can be performed with the aid of the replica method. For systems not containing 
quenched disorder one just needs to introduce replicas to average the logarithm in (Q). This consists in substituting 
log(Z) in (||) by Z m and evaluate (log(Z)) = lim m _>o((Z m ) — l)/m, where Z is the argument of the log in (||) and 
the angular brackets represent the average over the distribution of y. The formal procedure is similar to that used 
by Zippelius and coworkers to study vulcanization and the one used by Given and Stell to for liquids in random 
quenched matrix 13 1. The physical meanings of ours and their constructions is however very different. Both in the 
vulcanization and in the liquid cases the replica method is used to average over some kind of real quenched disorder, 
the random crosslinking occurring at the vulcanization transition in the first case and the quenched matrix in the 
second. In our case, there is no quenched disorder. The coupling with the reference configuration y is a theoretical 
tool that allows to probe regions of configuration space that have zero weight in the usual Boltzmann distribution, 
and that can allow us to give a description of freezing even in absence of quenched disorder. 

Writing explicitly the replicated partition function 



(Z m ) 



dy exp(- 



dxi...dx n exp 



a=l 



(4) 
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we see that the problem is reduced to an m + 1 component system in the limit m — > 0. As we will see in the next 
section, the procedure needs to be modified if the Hamiltonian H contains some quenched parameters, to take into 
account the denominator in (||). 

Let us now try understand qualitatively the behavior of Fq for small positive e and To < T < T c , in the simpler 
case where the two temperatures are equal. At e = the probability that the replica x would stay in a same local 
minimum of the replica y is exponentially small. While, when e > 0, the case in which the replica x stay near to the 
replica y is energetically favored. The system can therefore stay in two different phases 

• Replica x different from y (q very small) and its free-energy F(T, e) ~ F(T, 0). 

• Replica x near to y (here q w 1). The free-energy is given by F(T, e) » F(T, 0) — eq + TS(T). 

It is now clear that in this picture there is a first order phase transition at e w TS(T) with a discontinuity in the 
internal energy given by q. Moreover at e = the difference in free-energy among the two phases is exactly given by 
TS(T). For small e one finds that the transition line starts as T(e) = T c + Const, x e. 

The thermodynamic properties in the T — e plane (for different values of T') are quite interesting. The previous 
argument tell us something only in the region of small e, the fate of the first order transition for large e is a very 
interesting question. In principle such a computation could be done in structural glasses by using the replicated 
hypernetted chain approach of E3] and work is in progress in this direction p5| . As a first investigation we limit 
ourselves to study what happens in a generalized spin glass model, the spherical p-spin models with long range 
forces |fl6|| . We would like to stress another aspect that makes our approach interesting in connection with glass 
physics. Studying the usual Boltzmann measure of these models in the glassy phase one faces the problem that 
the configurations that give the dominant contribution at close but different temperatures look very different from 
each other (the so called chaotic temperature dependence of the measure). However, it happens in general that 
metastable states that dominate the measure at a given temperature T", that we call T'-states, remain metastable if 
the temperature is changed p7| . With our method, if we fix T" and we variate T we can "follow" metastable states 
in temperature In other words, we can modify the Boltzmann measure so as to give non-vanishing weight to the T'- 
states to different temperatures. This has important connections with cooling experiments in real systems. At a given 
cooling rate, the system equilibrates within the (super-cooled) liquid phase, until the glassy transition temperature 
T g is reached. This is the last temperature where the system is able to equilibrate. Below that temperature genuine 
off-equilibrium phenomena as aging and memory effects set in the system. However, one can expect that the system 
remains confined for a long time in the metastable state reached at T g . Indeed this has been observed in recent 
numerical experiments in Jl8[ |. For short enough times after reached the temperature T g the system will be found 
in local equilibrium in the "analytic continuation" of the state at T g . This hypothesis implies reversibility and can 
be valid only for times such that structural rearrangements can be neglected. In this perspective one would like to 
define restricted Boltzmann-Gibbs measures in which only the configurations with a given distance from the quenched 
configuration y have non-zero weight: 

P(x\y) = z^y) e ' mx)6iq{x ' v) ~ q) - (5) 

As the constraint on the value of q in (||) is a global one, the free-energy associated to the distribution (|5|), Vq(<?) = 
—T\ogZ((3,y), is related to Fq(c) by Legendre transform: 

Vq (q) = mm Fq (e) + eq. (6) 

This represents the minimal work required to keep the replica x at fixed overlap q from the replica y. In a situation 
with exponentially many minima, each one carrying vanishing contribution to the Boltzmann measure, one can expect 
a two minima structure of Vq at the mean field level. It is clear that a minimum should be found at the value of 
q characterizing the typical overlap among different states. This corresponds to having the second replica in one of 
the exponentially many equilibrium states, different from the one where the first replica lie. In addition, if one does 
not allow for configurations where q is spatially inhomogeneous there is a minimum corresponding to the two replicas 
globally in the same state.Q Therefore, for T — T' , the height of this secondary minimum with respect to the primary 
one has to be equal to the complexity S(T) multiplied by T. 



1 In real, finite dimensional systems the secondary minimum is washed out by the possibility of having configurations with 
inhomogeneous q but the same global overlap and smaller free energy then the homogeneous ones. This kind of configurations 
with "coexisting phases" lead to conclude that, as in usual systems, the potential has to be a convex function of q. 
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B. The annealed case 



A similar but different construction is the following: we consider two replicas with total Hamiltonian 

H e (x,y)=H(x)+H(y)-e £ w{ Xl - y k ) (7) 

i,k=l,N 

The difference with the quenched case is that for e ^ both x and y may not be equilibrium configurations. This 
construction, which we will see, gives similar results of the quenched one might be more advantageous in numerical 
simulations, where one can thermalize the two replicas at the same time. 

Let see what we do expect in this case if we make the approximation that q may be zero or 1. For T > T c there is 
a low q phase in which the free energy and the internal energy are just the same as at e = 0. The most interesting 
phase is the high q phase. At e = the free energy in the region T c < T < Td can be obtained for each of the two 
systems by minimizing the total free energy 

F(/)=2[/-T£(/,T)] (8) 

If we couple the two systems and we restrict the analysis to the pairs of configurations with q 1, we have that 

F(e,/) = 2/-TS(/,T)-e (9) 

In this case the minimum will be located at S(/, T) = also to temperatures higher than T c , i.e. up to a temperature 
T c 2 such that 



rp dE(J,T) 



dT 



(10) 

T c , 2 



In other word the presence of term proportional to e stabilizes the glass phase. As far as the difference in free energy 
of the liquid phase and of the glassy phase is of order (T — T c ) 2 , we expect that the second order phase transition is 
transmuted into a first order one with at a temperature T(e) = T c + Const, x e 1 / 2 and a discontinuity in the internal 
energy proportional to e 1 / 2 . This in contrast with the quenched case, where T(e) is linear in e. We see then, that 
although the symmetric coupling among replicas induces a kind of deformation of the landscape, so that e.g. the 
various transition temperatures are changed, the global situation is similar to that of the quenched case with the main 
difference in the behavior of the transition line. 



III. AN ANALYTIC COMPUTATION 



A. The p-spin spherical model 



The model that we are going to describe is the so called p— spin spherical model, which has become in the last years 
one of the simple reference mean-field models for the structural glass transition |l6| . It will be clear from the form of 
the Hamiltonian that this generalized spin glass is microscopically very different from a structural glass. The basis for 
its use as a model for the glass transition rely basically on the phenomenology of disordered models model with "one 
step replica symmetry breaking transition" , which is strongly reminiscent to that of structural glasses. This issue, 
that has been discussed widely in the literature will be the starting point for the application of our discussion 
to real glasses. Just to mention a few facts, the model presents an ideal glassy transition to a broken ergodicity 
phase with an extensive configurational entropy, to a low temperature zero complexity phase, very much like in the 
Gibbs-Di Marzio scenario, and its Langevin dynamics maps exactly in the schematic mode coupling theory ]l],Q] and 
its off-equilibrium generalization [fl9|-pl[ , capturing in this way many of the qualitative and some quantitative features 
of the supercooled and glassy relaxation. 
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Many study have shown explicitly that the random nature of the couplings in this models is not essential. In 
fact there have been found long range models with oscillatory couplings and the same basic physics of the model 



we are going to describe J^|23jjl9) . We than believe that many of the issues that we will discuss in this paper, and 
in particular the qualitatively features of the phase diagrams in the T — e plane are quite universal and reflect very 
general properties of the phase space. We can thus conjecture that the phase diagram for real glasses is similar to that 
of the generalized spin glasses if we only consider the order of the phase transition and the topology of the various 
transition lines in the T — e plane. 

The model is defined in terms of N real dynamical variables (spins) Si, (i — 1, ...,N) subjected to the constraint 
Si=i 3f = N and interacting via the Hamiltonian 



1,JV 



Hp— 2_, Jii,...,ipSii"S%p (11) 



ii<...<«. 



with independent centered Gaussian couplings Ji 1: ... : i v with variance Jf i — p\/{2N p v ). The model for p > 3 
has a one step replica braking, the transition pattern which has been repedetly shown to be deeply related to the 
Gibbs-Di Marzio entropy crisis mechanism PJ2^|. 

In spin models the natural way to couple two replicas consist in adding to the Hamiltonian a term — SiS'i an d 
with the usual definition of the overlap q — N^ 1 J2i=i n The overlap q is equal to one if the configurations of 

the two systems coincide. 



B. The Quenched case 

The quenched potential for the model is obtained inserting the Hamiltonian ( ]TT| ) in the the general definition ^ 
The two replicas potential is 



„ (TT , > JdS'eM-P'H(S'))F(T,e,S') 



JdS'exp(-(3'H(S')) 
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F(T,e,y)^(NP)- 1 \n\ dSeM-pH(S)+(3e £ S k S' k } 



i,k=l,N 



where the bar denotes the average over the J's. In addition to the m replicas needed to average the logarithm 
log(Z) = \im m ^ {Z m — l)/n, in order to average over the quenched parameters we also need to represent the 

denominator 1/z in (|l2] ) as lim„_ > o z n ~ x . As usual the computation is performed continuing from integer n and to. 
We have then n unconstrained replicas Sf , (a = 1, n) and to constrained replicas Sf, (a = 1, ...,to). The average 
over the quenched disorder induces a coupling among replicas, and as usual in this models, the order parameter for the 
theory is the matrix of the overlaps among all the n + m replicas of the system. This can be conveniently arranged in 
three matrices describing respectively the overlaps of the unconstrained replicas, Q a ,b — 1/N Y^i SfS!?, (a, b = 1, n) 
the overlap among the constrained replicas R a .p = ^/N^ i S"S^ , (a, (3 = 1, ...,m) and the mutual overlap among 
constrained and unconstrained replicas P a , a — V-^Si SfSf . In terms of these order parameters one finds 



(n,m) 



1,771 



a , b 



a,/3 

™ I /DP 

+2/3e^P lQ + -Trlog( ^ T ^ 



(12) 



where we have written f(q) = q p /2. From (|T^) one has to extract the terms proportional to to in order to evaluate 
the effective potential. The procedure that we have sketched here has been discussed in ref. p, p5| , |l7|| , to which we 
address the interested reader for the details. In the following we will concentrate to temperatures T > T s . In this 
region the matrix describing the unconstrained replicas has the simple form Q a ,b — 5a,b, which is valid both in the 
paramagnetic phase for T > T c and in the non ergodic phase for T s < T < T c with exponentially many states with 
vanishing weights. 




0.55 0.6 0.65 0.7 0.75 



FIG. 2. Phase diagram in the e — T plane for p — 4 and T 1 = T. The upper curve is the spinodal of the low q phase, the 
lower one the spinodal of the high q state, and the middle curve the coexistence line. The coexistence line touches the axes 
e = at T = T s , while the spinodal of the high q phase touches it at T = T c . For T < T s the spinodal of the low q phase 
remains finite and touches the T = axes at finite e. 
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In this conditions the constrained replicas have non zero overlap only with replica a = 1 and therefore, P a . a = 5a,i1- 
For the matrix R the most general ansatz needed here is a "one step broken" structure Jl0|]. Both the structure 
and the physical meaning of this ansatz have been discussed widely in the literature, and in standard notations we 
parameterize it by the three parameters qo,qi, x. 

In the following we will study the phase diagram of the model in the e — T plane in two situations: a) T' = T, 
corresponding to restricting the partition sum to the vicinity of a particular equilibrium state at each temperature, 
b) T" fixed, corresponding to probe the evolution of the free-energy landscape in the vicinity of a fixed configuration 
of equilibrium at T when T is changed. The Legendre transform of Fq(T, T',e), V Q {q,T,T') = min e Fq (T, T',e) + eq, 
admits the following expression in terms of the variational parameters defined above. 



V Q (q) = ~ |2/J/37(«) - P 2 ((1 - x)f( qi ) + xf(q )) + ^ln(l - Ql) 



-ln(l - (1 - x)q x - xqo) 
x 



qo 



1 - (1 - x)qi - xq 



(13) 



where Vq has to be maximized with respect to qo, q\ and x. Depending on the values of /?,/?' and q, the solution of 
the saddle point equations can be either replica symmetric with x — or x — 1, or display replica symmetry breaking 
with q\ ^ qo, and x + 0, 1 [ jl7| . 

We see from fig. [j] that the shape of the function V is the characteristic one of a mean-field system undergoing a 
first order phase transition. At high enough temperature Vq is an increasing and convex function of q with a single 
minimum for q = 0. Decreasing the temperature, we find a value T cr , where for the first time a point q cr with 
Vg(g cr ) = appears. The potential looses the convexity property and for T < T cr a phase transition can be induced 
by a field. A secondary minimum develops at T c , the temperature of dynamical transition signaling the presence of 
long-life metastable states. The minimum of the potential has received a dynamical interpretation in |5|,^6|,[l^] where 
it has been shown that its characteristics (internal energy, self-overlap, etc.) correspond to the states reached at long 
times by the evolution at temperature T starting at an initial time from an equilibrium configuration at temperature 
T' . In figure [l] we show the shape of the potential in the various regions. 




0.55 0.6 0.65 0.7 0.75 

FIG. 3. Latent heat of the transition as a function of the temperature for p — 4 and T 1 — T. 



Although the behavior of the potential function is analogous to the one found in ordinary systems undergoing a 
first order phase transition the interpretation is here radically different. While in ordinary cases different minima 
represent qualitatively different thermo dynamical states (e.g. gas and liquid), this is not the case in the potential 
discussed here. In our problem the local minimum appears when ergodicity breaks, and the configuration space splits 
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into an exponentially large number of components. The two minima are different manifestations of the existence of 
exponentially many states with similar characteristics. In the high q minimum the system S is in the same microscopic 
state as S', while in the low q minimum it can be in any other state among the exponentially many that contribute to 
the Boltzmann measure. The height of the secondary minimum, relative to the one at q = measures the free-energy 
loss to keep the system in the same component of the quenched one. At equal temperatures T = T' this is just equal 
to the complexity £ multiplied by T and it goes to zero at T s , where coexistence in zero coupling takes place without 
a release of latent heat and both minima lie at the same height. For T ^ T" the height of the secondary minimum 
also takes into account the free-energy variation of the equilibrium state at temperature T" when "followed" (i.e. 
adiabatically cooled or heated) at temperature T. 

The presence of the field e adds finite stability to the metastable states, and the transition is displaced to higher 
temperatures. The position of the transition line can be computed via the Maxwell construction. In figure we 
display the phase diagram of the p = 4 model in the case T" = T. The coexistence line departs from the axes e = at 
the transition temperature T s and reaches monotonically a critical point (T cr ,e cr ). For the mean field model under 
study one can see that the exponents characterizing the critical point are the classical ones. In figure || we also show 
the spinodal of the high q solution, which touches the e = axes at the dynamical temperature T c , and the spinodal 
of the low q solution for temperatures larger then T s . While the transition in zero field is not accompanied by heat 
release, a latent heat is present in non zero e. In figure || we show, in the same conditions of fig. 0, the latent heat 
Q = E+ — E- — e(q+ — q~) where E+ (q+) and E— (<?_) are the internal energies (overlaps) respectively of the high and 
low q phases. Notice that the latent heat vanishes at the critical point (as it should), and at T = 0. The coexistence 
line for T' fixed, in the interval T c < T' < T s is qualitatively similar to the one of figure |^ at high enough temperature, 
but (for T" > T s ) it never touches the axes e = 0. 




FIG. 4. Phase diagram in the e — T plane for p — 4 and T' — 0.523. The upper curve is the spinodal of the low q phase, 
the lower one the spinodal of the high q phase, and the middle curve the coexistence line. 

Even at zero temperature there is a first order phase transition in e, reflecting the fact that the ground state of the 
system is lower then the energy of the reference state when followed at T — 0. This can be seen in figure [| where we 
show the phase diagram for a value of T' such that T c < T' < T s . In figure || we see that the latent heat for T' =/= T is 
qualitatively similar to the one for T" = T at high temperature, while at low temperature it change sign and becomes 
zero only at T = 0. 
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FIG. 5. Latent heat of the transition as a function of the temperature for p and T' as in fig. ^. 

The high q phase roughly reflects the properties of the equilibrium states at temperature T' followed at temperature 
T, while the low q phase reflects the properties of the true equilibrium states at temperature T. At high temperature 
the high q phase is energetically favored, while at low temperature it has an energy higher then the one of equilibrium. 
The point where Q changes sign reflects this fact, and does not correspond to a second order phase transition. Finally, 
in figure]^ we show, for a fixed temperature the curve of q{e) obtained by the Maxwell construction. 

In closing this section we would like to comment on the use of the Maxwell construction for finite dimensional 
systems. In ordinary systems, the justification of the Maxwell construction is in the phenomenon of phase coexistence. 
Here we do not know what a supposed coexistence in physical space of the high q and the low q phase would mean. 
This, together with the related problem of of finding solutions with inhomogeneous q to the mean-field equations in 
finite dimensional models is an open problem and we let it for further investigation. For the time being we limit 
ourselves to look for support for our construction in numerical simulations. This will be the aim of next section. 









0.2 0.4 0.6 0.8 1 



FIG. 6. Equation of state for p = 4, T' = 0.523 and T = 0.609. The horizontal line corresponds to coexistence and is 
obtained by the Maxwell construction. 



C. The Annealed case 



The aim of this section is to show rapidly that the annealed potential gives results qualitatively similar to the 
quenched one. The detailed analysis can be found in |2^J. In this case we want to compute the average of the free 
energy of the two coupled replicas over the quenched couplings. The procedure to do this with replicas is explained in 
p8| , p7| . In this cases there are the overlaps describing the system can be conveniently arranged into 3 replica matrices 
Q a b = l/N^SfS^ Q' a b = l/N^S^Sf and P a ,b = All the matrices have the dimension n x n 

and the symmetry among the two system implies Q — Q' and P a j, — Pb, a - The most general ansatz needed in the 
problem is a "one step replica symmetry breaking" form where the matrices Q and P are parameterized respectively 
by the parameters (<?i,go = 0,x) and (p,pi,po — 0,x). The resulting free-energy is 

V A (q) = ~(3 [/(l) + f(q) -(1-x) (f( Pl ) + f( Ql ))} 

~ ~) [togC 1 +Pi - 9 - ?i ) + 1 °g( 1 - Pi + Q- 3i)] 



1 



[log(l +pi -q-qi + {-pi + qi) x)+ 



2(3x 

log(l -pi+q-qi + (pi+qi) x) 



(14) 



which has to be optimized with respect to qi,pi and x. In figure]^ we see that qualitatively the situation resembles 
to the quenched case, with a non convex potential function at low enough temperature. It is clear that the essential 
features of the phase diagram or the previous section with a first order transition line terminating in a critical point 
are present also in this case, although the actual values of the different characteristic temperatures will be different. 
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The annealed potential Va(q) as a function of q for p — 4 and various values of T, in order of decreasing temperatures 

from top to bottom. For p = 4 one has T c — 0.503 and T s = 0.544. The curves are normalized to have Va(0) = 0. 





FIG. 7. 



IV. REAL GLASSES 



A. The Model 



We have tested the ideas presented in this note for binary fluids. Although we did not perform a systematic study 
of the phase transition line, we found evidence for a first order transition in presence of coupling among replicas both 
in the annealed and in the quenched case. The model we consider is the following. We take a mixture of soft particle 
of different sizes. Half of the particles are of type A, half of type B and the interaction among the particle is given by 



H(x) = £ 

i<k 



a(i)+a{k) 



Xk\ 



12 



(15) 



where the radii (c(i)) depend on the type of particles. This model has been carefully studied studied in the past M. 
It is known that a choice of the radii such that obI&a — 1-2 strongly inhibits crystallization and the systems goes 
into a glassy phase when it is cooled. Using the same conventions of the previous investigators we consider particles 
of average diameter 1. More precisely we set 



a\ + 2{a A + (TBf + a\ 



= 1. 



(16) 



Due to the simple scaling behavior of the potential, the thermodynamic quantities depend only on the quantity T / p, 
T and p being respectively the temperature and the density. For definiteness we have taken p = 1. 

This is one of the simplest models of glass forming materials and we have chosen it because of its simplicity. The 
model as been widely studied especially for this choice of the parameters. It is usual to introduce the quantity T = /3 4 . 
The glass transition is known to happen around T = 1.45 Q|. It has been shown that aging appears below this 
temperature, and that the the time dependent correlation and response functions are well in agreement with the 
prediction of one step replica symmetry breaking below this temperature || . 

Our simulation are done using a Monte Carlo algorithm, which is easier to deal with than molecular dynamics. 
Each particle is shifted by a random amount at each step, and the size of the shift is fixed by the condition that the 
average acceptance rate of the proposal change is about .5. Particles are placed in a cubic box with periodic boundary 
conditions. 

Following the discussion in the introduction, we have introduced two copies x and y of the same system e have 
introduced the quantity q defined as 
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q(x,y) = j^^2 w ( x i -Vk), (17) 

where the sum over i and k runs over all possible N 2 /2 pairs of particles of the same kind. The function w is chosen 
in such a way that the quantity q counts the percentage of particles such that there is a similar particle nearby in the 
other configuration. The form we have considered is 

a 12 

W ( X ) = X 12 + fl 12 ' ( 18 ) 

with a = .3 The function w is very small when x » .3 and near to 1 for x < .3. The value of q will thus be a number 
very near to 1 for similar configurations (in which the particles have moved of less than a) and it will be much smaller 
value (less than .1) for unrelated configurations. The value of a has been chosen in such a way that q is insensitive to 
thermal fluctuations. 



B. Numerical results for the quenched case 

In this case the replica y is at equilibrium with the Hamiltonian H(y) while the Hamiltonian of the replica x is 

H(x\y) =H(x)-Niq(x,y) (19) 

where (3e = e. 

Our aim would be to find out if there is a first order transition in the plane e — T and to locate the transition line. 
In principle it is rather difficult to find out the precise position of a first order phase transition. The reason is quite 
simple: in a dynamical simulation the mean life in a metastable phase is exponentially large just near the transition. 

We have first followed an exploratory approach by monitoring the properties of the system as function of the 
coupling, temperature and their time variation, when we go from one phase to the other. In this case the results will 
be function of the speed at which we change the parameters and the precise value of the phase transition point will 
be obtained only in the limit of zero speed. 

This gives an approximate information on the position of the transition. We have taken a system with 66 particles 
and we have thermalized it at a given value of T for 2 k Monte Carlo sweeps (we have data for k = 7 — 17 in order to 
estimate the k dependence). 




0.2 0.4 0.6 0.8 1 1.2 



FIG. 8. The overlap as a function of the coupling for various values of F, from bottom to top F = 1.25, 1.30, 1.35, 1.40, 1.45. 
The number of particles is 66, the value of q is averaged over 10 different simulations, and the number of Monte Carlo step for 
each value of e in each simulation is 2 17 . 
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At this point we start the evolution of a second system which is coupled to the configuration reached by the first 
one (y) as in with a value of e = 1.2, taking as initial condition the quenched configuration y. We let then the 
system thermalize for other 2 k steps and we measure the overlap in the last quarter of the run. Starting from the 
final configuration we decrease the value of e by 0.1 and we perform the 2 k Monte Carlo sweeps. This procedure is 
repeated up to e = 0. 

In figure || we present the data relative of this procedure for k = 17. We see that at high temperature (low T) 
the value of q drops to nearly zero at already high values of e, while for lower temperature, it persists to high values 
down to low e. Data at lower value of k show a much smother dependence on e. The data are compatible with the 
possibility that for larger systems and for longer thermalization time a real discontinuity develops. This point deserves 
to be analysed in much greater detail. On the basis of the date of figure ||, in figure g we give a rough estimate the 
transition line in the plane T — e as the line where q = 0.7. 

There are two alternative methods which should be give more accurate estimates of the phase diagram. 

• We start from mixed initial conditions, i.e. half of the system in the phase with high overlap, half of the system 
in the phase with low overlap and we study which of the two phases becomes asymptotically stable. 

• We compute the free energy in each of the two phases (apart from a constant) by computing the internal energy 
and the overlap along a path that start form a fixed reference point in the e — T plane up to the final point. 




1.24 1.26 1.28 1.3 1.32 1.34 1.36 1.38 1.4 1.42 1.44 1.46 



FIG. 9. Rough estimate of the transition line in the e 
the previous figure is equal to 0.7. 



r plane, obtained as the value of e at which the overlap shown in 



In this note we explore only the first procedure. We fist thermalize the system for 2 k iterations. At this point we 
start the evolution of a second system which is coupled to the configuration reached by the first one (y) as in jl9| ) 
with a value of e = 4, taking as initial condition the quenched configuration y. The only difference form the previous 
case is that we substitute in the Hamiltonian (|l9|) q(x,y) by 



Qaix.y) 



Vk)G{x\) 



(20) 



where a^is the first component of the Xi and G{x) = —1 for x < L/2 and G(x) = 1 for x > L/2. In this way we force 
the overlap of the particles in the first and second half of the box to small and large values respectively. In this way 
we have prepared the starting point of the runs done (for others 2 fe sweeps) at different values of e with the usual 
Hamiltonian (119). In figure hO we present the data for a system of 512 particles for T — 1.35. We see clearly that for 
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high values of e the system evolves towards high values of q while for small e, q decreases to low values. The value of 
e that separates the two situations can be estimated to be around e = 0.9 in this figure. 

The dynamics of equilibration from that initial conditions presents interesting features. In figure [□] we show the 



overlap profile at times tk = 2 fe for e = 1.5, a value such that the high q phase is stable. Figure 12 shows the same 
thing in a case, e = 0.3, where the low q phase is stable. We notice that differently from the usual cases of first order 
transitions where the dominating phase grows at the expenses of the metastable one via a surface mechanism, here 
the dominant dynamics seems to occur in the bulk. The investigation of this kind of dynamics certainly deserves more 
attention than the one that we have dedicated to it here. 




10000 



100000 



1e+06 



FIG. 10. The value of the overlap for different values of e (from top to bottom e — 1.5, 1.2, 0.9, 0.6, 0.3), as a function of 
Monte Carlo time in a logarithmic scale. 




14 



FIG. 11. Time evolution of the space dependent overlap for a system prepared in the low overlap phase for < x < 4 and 
in the high overlap phase for 4 < x < 8. The number of particles is 512, and the box size is 8. The different lines represent the 
density profile averaged over y and z at different times tk = 2 k , for k = 0, 17. The values of F and the coupling, respectively 
F = 1.35 and e = 1.5 are such that the stable phase is the one with high overlap. The curves corresponding to higher times are 
higher in the low x region. 




FIG. 12. The same as in figure [ll] but with a value of the coupling (e = 0.3) such that the stable phase is the one with low 
overlap. The curves corresponding to higher times are here lower in the high x region. 



C. Numerical results for the annealed case 



In this case we have two replicas that evolve in parallel with Hamiltonian 

H(x, y) = H(x) + H(y) - Neq(x, y). 



(21) 



To support our first order transition picture we present hysteresis data during temperature cycles. The procedure 
consists in first cooling a system in which the two replicas start form independent random condition at a low value of 
r, when a maximum value of Y is reached, the two configurations are set equal (and equal to one of the two) and the 
temperature is raised again. In figures [l3| and |lj we show data corresponding to e = 0.2 and e = 0.8 respectively. We 
see in both cases that the low q phase seems to be metastable for all the probed values of the temperature. On heating 
the system passes from the high q phase to the low q phase with a sharp discontinuity. In figure [IB] we present data 
for the internal energy for the same cycle of figure [l4|. The hysteretic behavior found there is an important indication 
of a first order phase transition. Notice that the high q phase, where the entropy is lower, has a lower internal energy 
and is present also in the liquid phase for T < T c . Indeed the difference in energy starts to be present just around 
T c . This behavior is in agreement with the theory. Low energy metastable states exist also for T > T c , but they have 
zero Boltzmann weight in that region. Coupling two replica together enhance their probability so that their existence 
can be observed [f|i]||||l|,|§ . 
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FIG. 13. Overlap q for two coupled replicas in the annealed case, during a temperature cycle at small coupling e = 0.2 as a 
function of V. The cooling rate is 2 15 step for value of F, the number of particle is 66. We start from high temperature F = 1.3 
we cool down to T = 1.6 (lower curve). At this temperature we set the two configurations equal to one of them and we heat 
again (higher curve). 
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FIG. 14. Overlap q during a temperature cycle as in the previous figure, at higher coupling e = 0.8 as a function of T in 
the region F = 1 — 1.6. 
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FIG. 15. Internal energy E during the same temperature cycle as in the previous figure. The high curve corresponds to 
cooling, the lower one to heating. 

V. CONCLUSIONS 

In this paper we have given theoretical arguments and numerical support in favor of a description of the ideal glassy 
transition as the limiting point of a first order transition line in the plane of the temperature and the coupling among 
real replicas. 

Similar conclusions are obtained in the quenched construction and the annealed one. The first construction studies 
the implications of a Boltzmann-Gibbs distribution limited to the configurations which have a fixed overlap with a 
quenched configuration. The second construction studies a Boltzmann-Gibbs distribution of two systems on the same 
foot and with fixed overlap. 

In both formalisms the glassy state is described as a state where there are two phases. A "confined phase" with high 
value of the overlap, and a "deconfined phase" where there is minimal correlation. We have shown that the glassy 
transition, which is Ehrenfest second-order for zero coupling, becomes first order as soon as a non zero coupling is 
introduced. A detailed computation in the example of the spherical p-spin model has shown that that there is a first 
order transition line in the e — T plane terminating in a critical point. This result should be robust beyond mean field. 
Indeed, we have discussed how the numerical simulations for binary soft-sphere mixtures at non zero values of the 
coupling e support the theoretical picture. A much greater numerical effort would however be needed to locate with 
precision the transition line and to study the interesting issue of the nature of the glassy critical point. A key point in 
our analysis, in going from infinite range to short range models, is the possibility to use the Maxwell construction to 
estimate the topology of the transition line. In ordinary systems the validity of the Maxwell construction is intimately 
related to the phenomenon of phase coexistence. In the case of glasses we do not know what a supposed coexistence of 
the confined and the deconfined phase would mean. Our numerical simulations, which seem to confirm the first order 
transition picture, indicate an overlap dynamics different from the usual domain growth in spinodal decomposition. 
This question, and the one of the decay of the high q phase in the metastable regime are deeply related to the problem 
of restoring of ergodicity when the barriers are finite. This is at the heart of the glassy physics and further numerical 
and theoretical effort is certainly needed. 
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